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A class of algorithms for detecting abnormally short-holding-time 
trunks has been developed that utilizes individual trunk data available 
in EADAS/ICUR {Engineering and Administrative Data Acquisition 
System/Individual Circuit Usage Recorder). This data consists of a 
two-dimensional statistic that compresses the raw trunk measure- 
ments — the state of the trunk (busy or idle) sampled every 100 or 200 
seconds — into a manageable form. Because this data is essentially a 
sufficient statistic for the stochastic process used to model the (unob- 
servable) trunk state measurements, one of the algorithms developed 
is Wald's sequential probability ratio test. Two of the algorithms de- 
veloped have been implemented in WAN (Individual Circuit Analysis 
Program), and are currently being used to test trunks associated with 
the No. 1 crossbar, No. 5 crossbar, crossbar tandem (1XB, 5XB, XBT), 
and step-by-step switching machines. The focus in this paper, however, 
is on the modeling and analysis aspects of the problem, and only slight 
attention is paid to the various trade-offs and real-world constraints 
encountered in implementing the algorithms. 



I. INTRODUCTION 

A message trunk, the basic connecting link in the switched telephone 
network, provides the communication path between switching machines 
as well as certain call setup capabilities, such as supervision, signaling, 
and ringing. For an important class of trunk faults that cause call failure, 
the trunk is released by the switching system upon customer abandon- 
ment and is again available to fail another call. As a result, a single un- 
detected faulty trunk of this type can fail a significant fraction of the 
offered attempts to a group and will characteristically have an abnor- 
mally short holding time. 

Because of their potential service impact, significant efforts have been 
made to understand and quantify the impact that such abnormally 
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short-holding-time trunks have on central office and network service. 1-4 
It is now widely understood as a result of these studies that this generic 
trunk fault results in a fraction of service attempts "killed," which is out 
of all proportion to their number in the trunk population. Consequently, 
traffic data available from new and existing traffic data-acquisition 
systems has been viewed in the light of increasing interest in trunk-fault 
detection. In particular, with the advent of the Bell System EADAS/ICUR 
(Engineering and Administrative Data Acquistion System/Individual 
Circuit Usage Recorder) system, 5 it was natural to ask whether the new 
individual trunk data available could be used to detect such "killer" 
trunks.* 

This paper discusses the theoretical aspects of a class of killer-trunk 
detection algorithms that utilize the individual trunk traffic data 
available in EADAS/ICUR. These algorithms were designed for, and 
practical versions of them are presently implemented in, the ICAN+ 
portion of the EADAS/ICUR system. We focus here, however, on the 
problem formulation, modeling, and analysis aspects of the algorithms 
without bringing in many of the diverse factors and trade-offs encoun- 
tered in the actual implementation. 

Because the holding time of a trunk affects the statistical properties 
of the trunk data in EADAS/ICUR, it is natural to formulate the killer- 
trunk detection problem as a problem in the testing of statistical hy- 
potheses. In this context our modeling effort is basically an attempt to 
precisely define the state of a trunk (normal or killer) and expose the 
relevant underlying distributions. Well-known aspects of the theory of 
hypothesis testing are then applicable and immediately suggest a number 
of different tests. Sequential tests are naturally considered since the 
EADAS/ICUR data evolve sequentially in time. Questions about the ro- 
bustness of the models, and the structure and performance of statistical 
tests, are addressed using standard analytic tools. 

The material in this paper has been organized into six major sections, 
whose content we briefly describe. After considering the data available 
in EADAS/ICUR (Section II), we proceed to model a trunk (Section 3.1), 
motivate an appropriate set of statistical hypotheses suitable to our 
problem (Section 3.2), and briefly review several classical tests for de- 
ciding between statistical hypotheses (Section 3.3). With these prelim- 
inaries out of the way, we develop individual trunk algorithms based 
solely on individual trunk data. Proceeding in a heuristic manner, we 
use the individual trunk data to "derive" an ad hoc killer-trunk-detection 



* The term killer trunk has been widely adopted in referring to a faulty switching- 
machine-accessible trunk in a group whose average holding time is substantially smaller 
than the average group holding time. 

f Individual circuit analysis program — a software program that analyzes much of the 
EADAS/ICUR traffic data and maintains the system data base. 
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algorithm (Section 4.1). Although the insight gained in proceeding in 
a heuristic fashion is significant, we shift our emphasis in Section 4.2 and 
rigorously derive an optimal test statistic. It is interesting to find that 
the ad hoc statistic is essentially one of two symmetric statistics which 
comprises the optimal test statistic. The relationship between these 
individual trunk statistics is further explored in Section 4.3. 

In Section V we factor grouping information (which essentially 
identifies all trunks common to a trunk group) into the picture, and 
develop detection algorithms tailored to trunks associated with the 
No. 5 crossbar switching machine. This development necessitates 
modeling the 5XB trunk-group selection procedure, and several results 
due to Forys and Messerli 2 are utilized. In Section VI we shift our dis- 
cussion to the performance of the 5XB group algorithms, deriving ap- 
proximate expressions for the mean statistic update and mean detection 
time in Sections 6.1 and 6.2, respectively. The paper concludes in Section 
6.3 with an approximate analysis of the false-alarm probability of the 
5XB group algorithms. 

II. EADAS/ICUR DATA 

The structure of a killer-trunk detection algorithm is largely de- 
pendent on the type of individual trunk measurements available.* In 
EADAS/ICUR, the raw (unobservable) data consists of the state of each 
trunk (busy or idle) every 100 or 200 seconds. Fortunately, the data ac- 
cumulations available essentially summarize all the relevent information 
in the raw data. 

2. 1 Switch-count and transition data 

The data available from the EADAS/ICUR system, which can be used 
to distinguish between normal and killer trunks, is obtained by sampling 
individual trunks at 100 or 200 second intervals. This data consists of 
periodic accumulations (typically hourly, two-hourly, or three-hourly) 
of both the Busy states, and the State transitions. The busy state accu- 
mulation is usually referred to as the switch count. For the 200-second 
sampling option with a one-hour accumulation period, the switch-count 
is an integer between and 18. A state transition occurs whenever the 
state of a trunk (busy or idle) is different at two successive scans. For the 
200-second sampling option with a one-hour accumulation period, the 
number of state transitions is an integer between and 17. 

If we denote the ith scan during an accumulation period in which m 
scans occur by x,, and let and 1 correspond to trunk idle and trunk 



* Until very recently, almost all trunk-traffic measurements were obtained on a group 
rather than on an individual trunk basis. 
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busy, respectively, the available data may be written 

m 

(i) n{m) = £ xi (switch count) 

m 

(ii) t(m) = Yl \*i — *i-i| (state transitions). 

i = 2 

Thus, the raw (unavailable) data in the form of a binary sequence, 

x m — \ x Ij' ' '> x m)> 

is compressed into the two statistics n(m) and t(m). 

Because the holding time of a killer trunk is generally on the order 
of a few tens of seconds, it should have substantially more state transi- 
tions than a normal trunk, for a given switch count. Figure 1 illustrates 
the sampling process on both a normal and killer trunk.* For the pur- 
poses of this figure individual calls are represented by rectangles, call 
durations correspond to the width of the rectangles, and a half-hour 
accumulation period with the 200-second sampling option is used. 

We note in passing that for the 200-second sampling option, very little 
information is lost by "compressing" the raw data x m = (x-y,- • ;x m ) into 
the two statistics n(m) and t(m). Thus, normal conversation lengths tend 
to be in the vicinity of 3 to 4 minutes and, hence, with the 200-second 
sampling option, we expect that only adjacent samples are significantly 



* The realizations shown in Fig. 1 are more or less typical for a 5XB trunk group with a 
mean group holding time of approximately 4 minutes operating at about 40-percent oc- 
cupancy, and having a killer trunk with a mean holding time of approximately 1 min- 
ute. 
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correlated. If |*,}™i is Markovian, then [n{m),t{m)] is almost a sufficient 
statistic ([n(m),t(m),x\,x m ] is sufficient) for x m (see Section 4.2). Note 
also that for a trunk in the killer state, succesive samples should be es- 
sentially independent (for both sampling options). 

2.2 Grouping Information 

In addition to the switch-count and transition data available from the 
EADAS/ICUR system, we are also able to utilize a system map to identify 
(i) all trunks common to a trunk group, and (ii) the trunk-selection 
procedure* associated with the trunk group. It turns out that using this 
grouping information, 1 in addition to the switch-count and transition 
data, enhances the detection potential considerably. 

Thus, we divide the class of algorithms into two types according to 
whether or not grouping information is utilized. The first type, which 
uses only switch-count and transition data is referred to as an individual 
trunk algorithm. These individual trunk algorithms are applicable to 
all trunks — including two-way trunks — independent of the type of 
switching machine they are associated with. They do however assume 
knowledge of the trunks nominal holding time. The second type of al- 
gorithm uses the grouping information in addition to the switch-count 
and transition data and is referred to as a group algorithm. Group al- 
gorithms are "tailored" to a specific kind of trunk-selection procedure 
and, hence, apply to trunk groups associated with specific switching 
machines. For the purposes of this paper, the trunk-selection procedure 
considered is random selection of idle trunks. This procedure models 
the selection procedure of trunk groups associated with the 5XB 
switching machine. Group algorithms generally apply only to one-way 
trunk groups. 

III. PRELIMINARY CONSIDERATIONS 

In attempting to quantify the intuitive notion that a killer trunk will 
exhibit more transitions than a normal trunk (see Fig. 1), for a given 
switch count, it is natural to consider the transition probabilities: 

Pi.o = P|*i+t = 0/x t = 1J 
and 

Po.] = P\x t+T = Ux t = Oj, 



* The map in EADAS/ICUR indicates the type of switching machine that the trunk group 
is associated with, and this allows us to model the trunk-selection procedures (see 
Ref. 2). 

+ We will consistently use "grouping information" to refer to both the identification 
of all trunks common to the group and the trunk-selection procedure associated with the 
group. 
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where x t denotes the state of a trunk at epoch t and t denotes the sam- 
pling interval. Of course, to evaluate these transition probabilities, we 
must be concrete about how we model a trunk. 

Before tying ourselves down to any specific model, however, it is useful 
to view these conditional probabilities in a canonical form. Thus, suppose 
we begin by assuming only that the binary valued process x t is stationary. 
We have then the following simple result: 

Lemma 1 : Let x t be a binary valued, stationary random process and let 
p and /?(•) denote its mean and couariance function, respectively. 
Then, 

(i) Pi.o(p,r)=(l-p) ' 



1- 



m (1) 

(w) pPi,o(p,t) = (1 - p)P .i(p,t). (2) 

Proof: Part (i) is a consequence of the definition of i?(-). That is, 

R(T)=E(x t X t+T )-p 2 

= P(x t = l,x t+T =l)-p 2 , 

where 

p = E(x t ) = P(x t = I). 

Part («') follows from the two identities: 

P = pPi.i(p,t) + (1 - p)Po,i(p,t) 

and 

1-Pi,i(p>t)+Pi.o(p,t). 

A consequence of this result is that uncorrelatedness and independence 
are equivalent: 

Corollary 1 : For the process in lemma 1, x t ,x t+T are independent if and 
onlyifR(r) = 0. 

Note that the dependence of #(•) on p has been suppressed for conve- 
nience. 

3. 1 Modeling an Individual trunk 

A particularly simple way to model a trunk is as the server in a single 
server loss* system with a Poisson arrival process and an exponential 
service time distribution. This model is commonly denoted by 
M/M/1-loss. 6 Let x t denote the state of the server: 



* In a loss system, customers who are blocked depart without waiting. 
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1 if server is busy at epoch t 
if server is idle at epoch t 

and R(-) the covariance function of x t . It is easily shown 7 that for the 
M/M/1-loss system, 

R(T) = R(0)exp\-(X + n) T ), 

where X and p are the mean arrival and service rates, respectively. Thus, 
Pi,o(p> t ) m ay De written as: 

Pi,o(p,t) = (1 - P) jl - exp (f^ 1 ) ), (3) 

where the trunk occupancy p is equal to X/(X + p). Throughout this paper 
we will be concerned with t = 100 or 200 seconds and a nominal holding 
time 1/p in the vicinity of 3 minutes. The mean holding time of a killer 
trunk l//x* will always be expressed as 1/rp with r typically in the range 
5 to 15. Thus, if we denote \n by S, we may write the transition proba- 
bility P\x t + T = 0/x t = 1| for a trunk with mean occupancy p as 



Pi,o(p,r) = (1 ~ p) [l " exp (y^T")]. 



(4) 



where r = 1 corresponds to a normal trunk. (Since P\ t o ■ 1 — p implies 
that x T ,*2 T ,- ■ • are independent, we will assume independence for r 
sufficiently large in subsequent sections.) 

Figure 2 is a plot of Pi,o vs p corresponding to S = 10/9 (200-second 
sampling and a 3-minute mean holding time) for several values of r. P lf0 
is essentially equal to 1 - p for r > 5. Figure 3 is a similar plot of Pi i0 vs 
p corresponding to 100-second sampling and a 3-minute mean holding 
time (S = 5/9). In this figure Pi i0 is essentially equal to 1 — p for r ^ 
7.5. 

Before putting too much emphasis on the transition probabilities in 
Figs. 2 and 3, it is prudent to consider the effect of factoring more realistic 
assumptions into the single server loss model. Thus, while the Poisson 
arrival process assumption is probably a reasonable assumption for a 
trunk in a 5XB trunk group (random selection of idle servers), it poorly 
models the overflow nature of the traffic offered to trunks in a 1XB/XBT 
trunk group. + In the latter case, it is more appropriate to model the input 
stream to a trunk as a peaked process. 6 Figure 4 is a plot of Pi, vs p 
parameterized by the peakedness (z) of the input stream. This figure 
is based on an expression for P 10 derived for a GI/M/1-loss model* with 



t The trunk-selection procedure for 1XB and XBT trunk groups is essentially a two- 
sided ordered hunt. 2 

J GI/M/1-loss denotes a single server loss system with a renewal process input stream 
(GI) and an exponential (M) service time distribution. 
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Fig. 2 — 1 -*■ transition probability for the 200-second sampling option. 
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Fig. 3 — 1 — transition probability for the 100-second sampling option. 
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Fig. 4 — i -► o transition probability for the GI/M/1 -loss model with a switched Poisson 
arrival process — 200-second sampling option. 

a switched Poisson input stream (commonly used to model overflow 
traffic). 8 Appendix A contains several details on the model and deriva- 
tion. It is clear that the effect of peaked traffic on the transition proba- 
bility is very small (2 = 1 corresponds to a Poisson stream). 

Recent data 9 indicates that the service time distribution of a normal 
trunk has a coefficient of variation significantly greater than 1 (the ex- 
ponential case). Thus, in Appendix Awe derive the covariance function 
of the server process x t for an M/G/l-loss* model with a mixed exponential 
type of service distribution. Figure 5 is a plot of Pi, vsp parameterized 
by the coefficient of variation of the mixed exponential service distri- 
bution. We see that increasing the coefficient of variation has a noticeable 
effect on the transition probabilities, but the effect is to increase the 
discrimination between the normal and killer-trunk transition proba- 
bilities. 

Thus, it would appear that the transition probabilities based on the 
M/M/1-loss model are reasonably robust to perturbations in the trunk 
model. In addition, one suspects that using these transition probabilities 
in a detection scheme, which exploits the basic differences between killer 
and normal trunk transitions, might lead to a conservative design. 



1 M/G/l-loss denotes a single server loss system with a Poisson (M ) input stream and 
a general (G) service time distribution. 
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Fig. 5— 1 -* transition probability for the M/G/1 -loss model with a mixed exponential 
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3.2 Testing statistical hypotheses — a basic Idea 

Suppose a trunk has constant mean occupancy p (Ex t = p) and we 
observe it for h seconds during which n switch counts and tin 1 — *■ state 
transitions accumulate. We may write 



and 



where 



Hence, we have 



and 



n = L x kT 

k=i 



ho - Jl [Xkr ~ X(k-l)r] , 
k = 2 



h , f0 2^0 

m = — and z~ = I 

r ll 2<0 



E(n) = p fy 



(5a) 



E(t w ) = (m- l)E([x kT - x {k - l)T ]-) = (; - l) pPi.o- (5b) 
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Thus, for example, a trunk at 20-percent occupancy (sampled at the 
100-second rate) would accumulate 108 switch counts, on the average, 
in 15 hours. The mean number of 1 — state transitions in this time 
interval corresponding to a normal trunk is 43, but the corresponding 
mean number of 1 -> state transitions for a killer trunk (during the 
same time interval) is 86. (Referring to Fig. 3, we see that Pi,o(p = 0.2) 
is 0.40 and 0.80 for a normal and killer (r ^ 7.5) trunk, respectively.) 

With this example in mind, it is natural to consider a detection scheme 
of the following type: 

(i) Wait until we accumulate n switch counts on a trunk. 

(ii) When n = n , compare the accumulated 1 -»■ state transitions 
tio to some threshold Tf. 

(Hi) If tio > Tf, decide that the trunk is a killer, otherwise decide that 
the trunk is normal, [t l0 is not directly available (see Section 2.1), but 
it may be estimated by t/2.] 

If the trunk occupancy were known and fixed, this scheme would ap- 
pear to be very reasonable. The analogy to the usual scheme suggested 
for deciding between a fair and a biased coin is clear: no is the number 
of (hopefully) independent experiments (analogous to the number of 
coin tosses), with each experiment having just two possible outcomes: 
the scan which follows the switch count is either or 1. Thus, each switch 
count is associated with a 1 -* state transition ("heads") or a 1 -*■ 1 state 
transition ("tails"). 

From the point of view of statistical hypothesis testing, we are thinking 
of two underlying states: 

Null hypothesis H : Pi.o(p) = P trunk normal 
Alternate hypothesis Hi'. Pi,o(p) = P* trunk killer 

Thus, our intuition suggests that a threshold test of the type sketched 
above is natural for distinguishing between H and Hi. We will see 
(Section 4.1) that a (nonoptimal) test of this form arises naturally from 
pursuing the coin tossing analogy further. 

3.3 Problem formulation — sharpening the focus 

To simplify matters, assume to begin with that 

(i) The nominal mean holding time l/ix is known. 
(ii) The trunk occupancy p is known and constant. 
(Hi) The switch-count and transition data accumulations n(m) and 
t(m) are continuously available (scan by scan). 

With these assumptions, it is an easy matter to conceptually describe 
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and let 



the "optimum" scheme for deciding between the two simple hypothe- 
ses, 

Hq\ Trunk normal (mean holding time 1/**) 

if i: Trunk killer (mean holding time lAyi), 

with it understood that "trunk" refers to one of the specific models de- 
scribed in Section 3.1 (for concreteness assume the M/M/1-loss 
model). 

Let x m = (xi,« • ;x m ) be the sequence of trunk states up to and in- 
cluding the mth scan, and let the available data be (as before) 

m m 

n(m) = £ xi, t(m) = £ |x,- - x,-i|. 

Let 

Pim(n,t) = P(n(m) = n,t(m) = t/Hi) i = or 1 

Pom(n,t) 

The joint probability distributions P im (n,t), i = 0,1, are well defined 
for any specific trunk model, but they may be nontrivial to derive. £ m (-,-) 
viewed as a function of the vector random variable [n(m),t(m)] is re- 
ferred to as the likelihood ratio statistic and plays a central role in the 
theory of statistical hypothesis testing. More specifically, the optimum 
test (in a variety of senses) for deciding between two simple hypotheses 
involves suitably comparing £ m to a threshold (or thresholds) in order 
to make a decision. 

We briefly review two optimum tests, the Neyman-Pearson (fixed 
sample) test 10 and Wald's sequential probability ratio test (SPRT), using 
notation appropriate to our (discrete) problem. 

3.3. 1 The Neyman-Pearson test 

Suppose a and denote the type 1 and 2 errors* of the test, 

Choose Hi if £ m ^T 

Choose H if £ m < T, 

and suppose a' and $' denote the type 1 and 2 errors of any other test 
(requiring m samples) for deciding between H and H\. Neyman and 



* The type 1 and 2 errors, « and 0, are often referred to as the probability of false alarm 
and the probability of miss, respectively (a = probability of choosing Hi given Ho is the 
true state, = probability of choosing Ho, given Hi is the true state.) 
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Pearson's classical result is: if «' ^ a, then ft ^ 0. Thus, of all tests re- 
quiring m samples and having a false-alarm probability not exceeding 
a, the likelihood ratio test achieves the minimum probability of miss 
(maximum probability of detection). Since a = P(£ m > T/H ), choosing 
a sample size m and threshold T to achieve a ^ «o requires knowledge 
of the (conditional) distribution of £ m . Similarly, having chosen m and 
T, calculating = P(£ m < T/Hi) requires the distribution of £ m (con- 
ditioned on Hi). Note also that with such a fixed sample test, we decide 
in advance to accumulate exactly m samples before making a decision. 
In many contexts, data accumulates sequentially in time, and rigidly 
requiring m samples — independent of the particular realization that is 
unfolding — is not an optimal strategy. 

3.3.2 Wald's sequential probability ratio test 

Using Wald's SPRT, 10,11 we continue to update £k, k = 1,2,- • • and 
defer a decision as long as C k dTo,Ti). We make a decision the first time 
?k falls outside the interval (T ,Ti). Thus, 

ifTo<^<Ti k = 1,2,- • -,m - 1 

and e m <$ (T ,Ti), 

then choose H\ if £ m ^ T\ 

and choose H if £ m = T . 

Clearly the stopping time m of the SPRT is a random variable, and the 
mean of m (given either hypothesis) is a measure of the time it takes to 
reach a decision. (Under a wide variety of circumstances, the SPRT ter- 
minates with probability 1.) Let Ei(m) (i = or 1) denote the mean 
stopping time, given that hypothesis i is in effect. Given a SPRT with type 
1 and 2 errors a and /3, and with mean stopping times E (m) and -Ei(m), 
consider any other test (sequential or not) with type 1 and 2 errors a' 
and 0', and with mean stopping times E' (m) and E\(m). The SPRT has 
the following optimal character 10 

if ex' ^ a and & s 0, 

then E' (m) S E (m) and E\(m) ^ Ei(m). 

Thus a SPRT is superior to a fixed sample test, if both tests have the same 
type 1 and 2 errors, in the sense that on the average it reaches a decision 
more quickly (under either hypotheses). 

In sharp distinction to the fixed sample test, the thresholds Tq and 
T\ required to approximately achieve specified type 1 and 2 errors are 
trivially determined. 11 On the other hand, even determining the mean 
and variance of the stopping time is often a difficult chore. 
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In Section 4.2, we explicitly calculate the SPRT* for the simple hy- 
pothesis testing problem described at the beginning of this section. 
Before looking at this optimum test, however, we describe an ad hoc 
algorithm which is very robust and consequently attractive from a 
practical point of view. 

IV. INDIVIDUAL TRUNK ALGORITHMS 

A basic underlying assumption in this section is that the normal mean 
holding time of a trunk is known. Thus, if the algorithms in this section 
are designed relative to a normal mean holding time of 3 minutes, they 
will not discriminate between normal trunks having a mean holding time 
in the vicinity of 40 seconds,* and an actual killer trunk with the same 
mean holding time — both of these trunks will be detected as killer 
trunks. 

The rationale for studying this type of detection problem is two-fold: 
from the practical point of view the simplicity of implementation and 
general applicability* of these algorithms is attractive, and EADAS/ICUR 
can flag trunk groups which should not be studied by the killer trunk- 
detection algorithms (thus preventing false alarms on normal short- 
holding-time trunks). From the theoretical point of view, it was natural 
to consider this problem before factoring group information into the 
picture. 

Another modeling assumption used in this section (as well as in sub- 
sequent ones) is that the arrival process is stationary within data accu- 
mulation intervals, but the mean arrival rate may change arbitrarily from 
one accumulation period to another. Since we use equilibrium analysis 
(e.g., in calculating Pio) we assume, in effect, that equilibrium is achieved 
instantaneously. 

4.1 An ad hoc algorithm 

The essential idea of the test suggested in Section 3.2, is to decide on 
the state of a trunk by comparing the number of 1 -*■ state transitions 
(tio) to some threshold Tf, conditional on having accumulated a fixed 
number of switch-counts. We heuristically* proceed to derive such a test, 
using a standard likelihood ratio formulation, and explicitly take into 
account the time-variability of traffic. 

Let x m = (*!,• • -,x m ) correspond to the (unobservable) binary se- 



* Based on the M/M/1-loss model for a trunk. 

* Trunks in special-purpose trunk groups (credit checking, weather, etc.) will typically 
have mean holding times in the vicinity of 40 seconds. 

* The individual trunk algorithms can be used to test any trunk — regardless of the type 
of switching machine the trunk is associated with. 

* The distributional assumptions made in this section are intuitively motivated, but 
cannot be rigorously justified. We examine these assumptions carefully in Section 4.3. 
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quence of trunk states during an accumulation period in which m scans 
occur. Let t\o(m) and n(m) be the number of 1 -* state transitions and 
switch counts associated with x m . Denote the conditional probability 

P(tioim) = t/n(m) = n) 

for a normal and killer trunk by P(t/n) and P*(t/n), respectively. These 
conditional distributions depend, of course, on the trunk's occupancy 
and on the particular trunk model we have in mind. [P*(t/n) also de- 
pends on the killer parameter r.] However, for the purposes of the heu- 
ristic development of this section, we do not precisely define which trunk 
model we have in mind. 

Since each switch count is associated with either al-*0oral-»-l 
state transition with probabilities Pi i0 and P lfl = 1 — Pi i0 , respectively 
for a normal trunk, and since we expect successive transition events on 
a trunk to be essentially independent,* it seems reasonable to assume 
that P[tio(m) = t/n{m) = n] for a normal trunk is binomially distributed 
with parameters n and Pi,o. This same argument applies to a killer trunk. 
Denote the binomial distribution with parameters n and p by b(k;n,p) 
k = 0,' ■ -,n, where 



b(k;n,p) = Q p k (l - P) n ~ k . 



Thus, we may think of a trunk with occupancy p during an accumulation 
period as having a conditional distribution 

P[t 10 (m) = t/n(m) = n] = b[t;n,P lt0 (p,r)], (6) 

with r = 1 and r = r corresponding to the normal and killer states of the 
trunk. (Recall that Pi,o(p,r) is essentially independent of r for r ^r with 
r = 7.5 and 5.0 for 100- and 200-second sampling, respectively.) With 
these assumptions, we may think of testing the two simple hypothe- 
ses: 

H Q : P(t/n) = b(t;n,P h0 ) P 1>0 = Pi.o(p.l) 

Hi: P(t/n) = b(t;n,Pl ) P\ tQ = P h0 (p,r ). 

If the 1 — *• transition and switch-count accumulations for two successive 
and contiguous accumulation periods are (£i,ni) and (£2.^2) respectively, 
we assume that 

P(ti,t 2 /n h n 2 ) = P(ti/ni)P(ts/ns). 

The idea here is that the only dependence between the two successive 



* The idea is that if significant correlation extends only one or two scans back, then 
successive transition events (events "triggered" by switch counts) should be essentially 
independent. 
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bit streams x l m = (xi,- • -,x m ) and x£, = (x m+h - • ;X2m) is essentially due 
to the dependence between x m and x m +v 

Thus, if we denote the transition and switch-court accumulations for 
the j'th accumulation period (in which m, scans occur) by [tio(mi),n(rrii)] 
during which the trunk has occupancy p„ we have 

P[tio(mi) = ti,- • -.tioirrik) = t k \n(mi) = m,* • -,n(m*) = rik] 

= n b[tr,n i ,P li0 (pi,r)], (7) 

where pi i ; = 1,« • •,& are the occupancies for the k accumulation periods. 
If tk = (tw - ;tk) and n^ = («!,• • ;Uk) consider the likelihood ratio: 

€(tk/n k )- II— — J — — . (8) 

,=i b[ti;ni,Pi t o{Pi,l)] 

Denote the log likelihood ratio* log £(tk/nk) by i(tkfnk) and note 
that 

Htk/n k ) = £ S(tj/iij) f 
i=i 

where 

* /4 , v , b[t/;n t -,Pi, (p/,r )] 

^U/n,) = log— — — — . 

b[ti\ni,Pifi{pi,l)\ 

The expression i(ti/ni) can be written as 

HUM = a(pi)ti - a(pi)rii 



with 



and 



i \ i ! ~ Pi.o(p.I) ro v 

a(p) = log- — *- (9a) 

1 - Pi.oip.ro) 



«(p) = a(p) + log^^. (9b) 

Pi,o{p,l) 



Thus, we have 



Ht k /n k ) = E [«(p.-)t,- - a(p,-)n,-]. (9c) 

i'=i 

Unfortunately, the occupancy in the ith accumulation period (p,) is 
unknown and hence equation (9c) cannot be used as a test statistic. One 
obvious "fix" is to estimate p, by p, = n,/m,, where n, and m, are the 



1 1 > T iftgW) > g(T) if g is monotone increasing, so the tests/ > T and g(l) > g(T)are 
equivalent. 
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switch count and the number of scans, respectively, during the ith ac- 
cumulation period. In stationary traffic, the estimate 

t / » 
fa = E nj / L mj 

7=1 / ;=1 

would be used [<r(pi) = l/Vi a (fa) if my = m for all;]. 

Corresponding to the sequence of accumulations, (£i,n,-,m;) i = 1,2," • • 
we define r, and /?,-, i = 1,2,- • • by 

r, = a(fa)ti - a(fa)rii, fa = n./m,- (10a) 

and 

Ri = /?,_! + r, with fl = 0. (10b) 

Thus, we arrive at the sequential test: 

(i) Compute Ri, i '.— 1,2,* • • and defer making a decision as long as To 

< Ri < Ti. 

(ii) If i '. = k corresponds to the first accumulation period for which 

Ri $ (To.Tr), then 

Rk ^ T => Trunk normal 

Rk^Ti^ Trunk killer. 

If we ignore the fact that we are estimating p, by fa, and by assuming that 
the various assumptions made are valid (see Section 4.3), we identify the 
above test as Wald's SPRT and as such To and T\ can be calculated as 
follows: 11 to approximately achieve type 1 and 2 errors, a and 0, re- 
spectively, a + < 1, choose 

To = log (^) (10c) 

and 

7\ = log(^). (10d) 

Throughout this section, we have assumed that the 1 -* transitions 
(t\o) are available when, in fact, only the total transitions (t) are available. 
It should be clear that t io can differ from fc/2 by at most ±%. To be pre- 
cise, let t \o(m), toi(m) be the number of 1 -* and — »■ 1 state transitions 
corresponding to a bit stream x m = (*i, • • •, x m ). If n(m) is the switch 
count corresponding to x m , then we have 

n(m) = tio(m) + t n (m) + x m (11a) 

and 

n(m) = toi(m) + tu(m) + Xi, (lib) 
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2 



where t u (m) is the number of 1 -*■ 1 state transitions. Therefore 

iio(m) + x m = toi(m) + xi, 
which together with t(m) = tQi(m) + tioim) yields 

t 10 (m) = i t(m) + ( Xl ~ 2 Xm ) (12a) 

and 

toi(m) = ^t(m)-(^Y^y (12b) 

Thus, we can write the statistic update (eq. 10a) as 

[It is easy to show that E[{x\ — x m )a{p)] = 0.] 

We conclude this section with an interpretation of the statistic update. 
Rewriting the statistic update as 

r = (a- a)t\Q - a(n - ti ) 

and using eq. (11a), we obtain 

r= (a-a)ti -at u -ax m . (13) 

Now, from eqs. (9a) and (9b), it is clear that a > a > 0. Thus, each 1 -*■ 

transition is weighted positively (evidence of a killer) while each 1 -*- 

1 transition is weighted negatively (evidence of a normal trunk). This 
is an intuitive explanation of the fact that the random walk (eq. 
(10b)) 

Rk - Rk-i + fk 

has a positive drift if the trunk is a killer and a negative drift if the trunk 
is normal. 

The fact that the update assigns a negative weight (—a) whenever the 
last bit (x m ) is 1 uncovers a modeling deficiency. Recall that in eq. (6) 
we assumed 

P(t 10 (m) = t/n(m) = n) = b(t;n,P l0 ), 

even though x m = 1 can not contribute to an observable 1 — »• transition. 
In this way we effectively modeled in a bias towards making "trunk 
normal" decisions. We can easily correct eq. (6) by conditioning on 
whether x m = or 1, obtaining: 

P(tio(m) = t/n(m) = n) = (1 - p)b(t;n,P 10 ) + pb(t;n - 1,P 10 ). 

Now, proceeding as before in formulating the log likelihood ratio yields 
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a statistic R&, where 

and 

r k = ?k + 9*. 
where r* is defined by eq. (10a), 

! + /-A-Wi - A.V 1 \ 

. \l-Pk/\ n k /\l-P l .o(p k ,r Q )J 

q k = log : (14) 

1 + \l-fa)V ~ 'r7j (l-Pi.o(M,l)/ 

and fa, tk, and n^ are the trunks occupancy estimate, 1 — * state tran- 
sitions, and switch count, respectively, during the fcth accumulation 
period. 

Thus, we obtain our original test statistic with the correction term qk 
added on. Note that q k ^ 0, qk — as fa -*■ 0, and qk — a as fa — 1. 
which is just the type of behavior expected, to offset the bias term in 

Having heuristically developed an ad hoc sequential algorithm that 
is intuitively appealing and easily implementable, it is natural to ask: 
how does it compare to the optimum sequential algorithm? In the fol- 
lowing section, we rigorously develop an optimum sequential test. 

4.2 An optimal algorithm 

Consider the two simple hypotheses: 
Ho: Trunk normal (mean occupancy p, mean holding time l/n) 
H\. Trunk killer (mean occupancy p, mean holding time 1 /r fi). 
The optimum test for deciding between the two hypotheses — in the sense 
of minimizing the mean decision time — for given type 1 and 2 errors, is 
Wald's SPRT (see Section 3.3), and it is based on the likelihood ratio 
statistic £ m (t,n) given by 

tmit,n)= — — — -. (15) 

P(t(m) = t,n{m) = n) 

Thus, it is clear that the ad hoc test described in Section 4.1 is not opti- 
mal, based as it is on an assumed conditional distribution, 

P(tw(m) = t/n(m) = n). 

Before proceeding to study eq. (15), we must define the trunk model 
precisely. In the developments that follow, we model a trunk as the server 
in an M/M/1-loss system (see Section 3.1). The model implies that the 
sequence of trunk states x t ,t = kr,k = 1,2,- ■ ■ is Markovian. Note that 
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although this appears to be a reasonable model for a normal trunk with 
200-second sampling, it ignores the conditional dependence "2 samples 
back," which is more important for 100-second sampling — e.g., x t given 
x t -r is independent of x t -2 T for the M/M/1-loss model. Taking this de- 
pendence into account in a trunk model would not be useful however, 
since the data needed to implement dependence "two scans back," is not 
available. 

Since we are modeling the sequence of trunk states as a binary valued 
Markov process x^, k = 1,2,- • •, in equilibrium,, ! it is clear that this process 
is characterized by = (Pi.o.Po.i), where Pi,o and Po.i are the transition 
probabilities 

P(x t+T = 0/xt = 1) andP(x t+r = l/x t = 0), 

respectively. (In general, a binary valued Markov process Xk T , k = 1,2,« • • 
in equilibrium, can be characterized by any two of the three quantities 
P, Pi t o,Po,i- For our special Markov process (based on the M/M/1-loss 
model), both Pi t o and Po.i and hence the process itself is determined by 
p alone.) Now having observed any m-tuple of the samples, which we 
denote by x m = {x\,- • %x m ), it is trivial to show that the statistic 

T(x m ) = (t(m),n(m),xi,x m ) 

is a sufficient statistic for 0. Thus, except for the initial and terminal 
states (x\ and x m ), the transition and switch-count accumulations 
summarize all the "relevant information" in x m . 

Our hypothesis-testing problem can now be formulated as follows: 
*i»*2,* • • is a binary- valued Markov chain in equilibrium with parameter 
= (Po.i.Pi.o) or 0* = (Po,i,P*,o)- That is, our two states are 

Hq: \xi] Markovian, characterized by 6 = (Po,i»Pi,o) 

H\\ \xi\ Markovian, characterized by 0* = (Po,i,Pi,o)- 

Now, because (t(m),n{m),x\,x m ) is a sufficient statistic for 0, we know 
that the likelihood-test statistic based on the raw (unobservable) data 
x m = (*i/ * m ,x m ) will be expressible in terms of t(m), n{m), x\ and x m 
only. Thus, instead of studying eq. (15), we proceed (for simplicity) to 
study the likelihood-ratio statistic: 

/(x m ) = log^^f . (16) 

P(*m) 

In Appendix B we study £ m (t,n) and find that it differs from £{x m ) only 
in an end-effect term. In $(x m ) this term depends on x\ and x m , whereas 
in i m (t,n) the corresponding term is a function of t and n. 
Since 
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P(x m ) = P( Xl ) U PiXi/Xi-,), 

/=2 



we may write 
so 



p(x m ) = p(xi)pa°pj»ps? i P{g o 

P*Ui) . . , Plo , , /l-Plc 



g(g w )"log p/ . +tiolog--+tiilog(- — — 1 
rUi) no \1 — i*io/ 

+ (o,log^l + toolog(7^). (17) 
Poi M — "01/ 

Note that a trunk with mean occupancy p is busy and idle with proba- 
bility p and 1 — p respectively, independently of the state it is in (normal 
or killer). Thus, 

P*(*i) 

log = and eq. (17) can be written 

PUi) 
i(x n ) = [(a - a)t 10 - at u ] + [(0 - b)t 01 - bt w ], (18) 

where the parameters b and are defined by 

1-Pr 



b = log (f^) (19a) 

\l — foi/ 

^ = 6 + log(^), (19b) 

and the parameters a and a are defined as in Section 4.1 (eqs. (9a) and 
(9b)). P,*)i and Pj i0 correspond to P ,i(p,r) and Pi,o(p,^) with r = r . 

Before discussing the symmetric structure of the optimum statistic 
[eq. (18)], we examine the P Uf i characteristics for the M/M/1-loss model. 
Using eqs. (2) and (3), we can obtain P ,i vs mean-trunk occupancy p for 
a normal (r = 1) and killer (r = r ) trunk. Figures 6 and 7 are plots for 
the 200- and 100-second sampling option, respectively, with a mean 
holding time of 180 seconds. It is clear from Fig. 6 that a -— 1 transition 
is just marginally more likely to occur on a killer trunk than on a normal 
trunk with a 200-second sampling rate. Although, the difference in the 
-*■ 1 transition probabilities between a normal and killer trunk in- 
creases substantially with the 100-second sampling rate, it is clear that 
these differences are still quite small — compared to the spread between 
the Pi,o and P\ t0 plots (see Figs. 2 and 3). Note that eqs. (2) and (4) show 
that 

1 - exp — — 

^ i = ^>lforr >l (19c) 

Po, l _ — s 



1 — exp 



1-P 
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Fig. 6— — 1 transition probability for the 200-second sampling option. 
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Fig. 7—0 — 1 transition probability for the 100-second sampling option. 
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and, hence, we have > b > 0. Using eqs. (4) and (19c), we see that 

P*o,i/Po,i - P\fi/P\,o and, therefore, 

- 6 = a- a (19d) 

Equation (18) shows that the optimum statistic is the sum of two 
symmetric statistics: 

(i) The statistic [(a — a)t\ Q — at\\\, which is essentially the ad hoc 
statistic (see eq. (13) and related discussion). 

(ii) An additional statistic [(0 — b)t i — fitfool which weights -»• 1 
transitions positively (evidence of a killer) and — * transitions nega- 
tively (evidence of a normal trunk). 

Note that by interchanging the role of and 1 in either of these two 
statistics, we obtain the other — b is obtained from a and is obtained 
from a by replacing Pi,o with Pq.i. 

By using eq. (11a) and the analogous equation 

m - n(m) = too(m) + t i(m) + x c m , {x% = 1 - x m ) (20) 

in eq. (18), the optimum statistic can be written 
Hx m ) = [atio(m) - an(m)] + [pt 01 (m) - b(m - n)] + ei(x m ), (21) 

where the end-effect term e\{x m ) is given by 

e\{x m ) = ax m + bx c m . 

To implement £(x m ) with only t(m) and n(m) available, necessitates 
estimating both ti Q (m) and t i(m) by t(m)/2. That is, using eqs. (12a) 
and (12b) in eq. (21) yields. 



/(x m ) = [«^-a n (m)] + {/3^-6[m- n (m)]] 

where 

e(xi,x m ) = (a - 0) ( 1 9 ) + e l(*m) 

*a-fr 



+ e(xi,x m ), 

(22a) 



= ( 2 ) (Xl + JCm) + 6t (22b) 



or 

>a + 



Hx m ) = ( CT „ ) £ ( m ) ~ (a ~ b)n(m) - bm + e(x lt x m ). 



(23) 



f Recall that a-p = a-b [eq. (19d)J. 
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In the development of the ad hoc algorithm, we assumed that the 
statistics corresponding to successive accumulation periods are inde- 
pendent. We conclude this section by examining the independence as- 
sumption and with some remarks on implementation. 

Thus, turning to multiple accumulation periods, suppose x' mi i = 
1,2,* • %k are the (unobservable) bit streams for k successive (and con- 
tiguous) accumulation periods, where x' m = (X(/-i)in+i»" • v£, m )- As- 
suming stationary traffic, and noting that \x, \f=i is Markovian, we can 
write 

Pixir ■ -,x* ) = n p(xi) x n i p( *; m+l/ *; m) (24) 



and, therefore, 
/UJm- • ;x k m ) = L i(x l m ) + Z log 

i'= 1 i= 1 



[P*(Xi m +i/Xi m )/P*(x irn+ i) 



(25) 



P\Xim+l/Xim )'P\Xim+l) 

where P(«) and P*(«) denote the distribution under Ho (trunk normal) 
and H\ (trunk killer), respectively. But, as we have seen, 

-ros 



P*(Xim+l/Xim) n 

= 1 — exp 



I-P 



(26) 



P*(x im+ i) 

and hence *; m +i and Xj m are essentially independent for ro sufficiently 
large. Therefore, 

P*(Xim+l/Xim) = P*(x im+ i) 

and, hence, eq. (25) may be written as 

k k—i 

i(x mt - • -,x k m ) = £ i(x' m ) - £ I{x im ;x im+1 ), 

i= 1 i'= 1 

where 

ii \ l *\Xim+l/Xim) /on\ 

I{x im ;x im+ {) = log— — — (27) 

"KXim+i) 

is recognized as the mutual information random variable, which plays 
a central role in information theory. 12 It is well known 12 that (under Ho) 
E\I(xi m ;xi m +i)\ is non-negative, and hence to ignore the end-effect 
term 

k-i 

H I(xi m ;xi m+ i) 
i=i 

by implementing the statistic 

i"=l 
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would tend to make a normal trunk look more like a killer trunk on the 
average. In Appendix C, however, we show that the mean end effect 
E{I{xi m ;x im +i)} is negligible compared to the mean statistic update 

The "optimal" sequential algorithm is implemented in the same 
manner as the ad hoc algorithm (see eqs. (10a) and (10b)) except that 
now, corresponding to the sequence of accumulations (t,,n,,m,)i = 
1,2,- • •, we define r, by 

'« + p\ . 

' l - — Orrii 



r ' = ( 2 ) tl ~ (a ~ b)m " ^ 



The term ti denotes the total number of transitions in the ith accumu- 
lation interval. As is the case for the ad hoc algorithm (which corresponds 
to b = ft = 0), the weights are functions of the trunk occupancy estimate 
Pi = ni/rrii. In a practical nonstationary environment, no claims of op- 
timally are made or implied. The term "optimal" is applicable only in 
the context of the equilibrium (e.g., stationary) model with known trunk 
occupancy. 

4.3 The ad hoc algorithm reviewed 

The assumption that tio(m) conditioned on the switch count n(m) 
is binomially distributed, is the basic assumption in the development 
of the ad hoc statistic. Although this assumption is incorrect (as we will 
soon see), the ad hoc statistic is essentially (except for an end-effect term) 
one of two symmetric statistics whose sum is the optimum statistic. Our 
purpose in this section is to examine the binomial assumption and to 
explain the relationship found between the ad hoc and optimal statis- 
tics. 

Since the optimal statistic was developed for a trunk modeled as a 
server in an M/M/1-loss system, it is natural to examine the binomial 
assumption (eq. (6)): 

P[tio(m) = t/n(m) = n] = fe[t;n,Pi.o(p,r)], (28) 

where Pi,o(p,r) is given by eq. (4) in this context. Consider a killer trunk 
with r sufficiently large and suppose x m = (*i/ • ;x m ) is the bit stream 
for a killer trunk during some accumulation period. Then, for all practical 
purposes [see eq. (26)], the trunk states xi i = 1,2,- • -,m are independent 
and identically distributed Bernoulli random variables: 

f P if x = 1 

P( Xi = *) = i - f n 

11 — P lfx = 0. 

Thus, it is clear that the switch-count distribution on a killer trunk is 
the binomial: 
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P(n(m) = n) = b(n;m,p). (29) 

Now suppose A m (t,n) denotes the number of binary m -tuples having 
exactly t 1 -*■ transitions and n ones. If x m = (xi,» • ',x m ) is a se- 
quence of trunk states for a killer trunk with n(m) = n, it is clear that 
each such sequence has probability 

P(x m ) = P n (l - p) m ~ n 

and therefore 

P(t 10 (m) = t,n(m) = n) = A m (t,n)p n (l - p) m ~ n (30) 

for a killer trunk (with r sufficiently large). Equations (29) and (30) show 
that for a killer trunk, 

(V)C) 

P(tw(m) = t/n(m) = n) , (31) 

C) 

where we have used the fact that 

A m (t,n).( m - n )Q. (32) 

It is interesting to note that while our assumed distribution for a killer 
trunk (28) differs from the correct distribution (31) — note that (31) is 
independent of p — there are some interesting similarities. For example, 
the assumed distribution peaks in the vicinity of (n 4- 1)(1 — p) and has 
mean equal to n(l — p) + whereas the true distribution peaks in the vi- 
cinity of (n + 1)[1 - (n/m)] and has mean equal to n[l - n/m]. Note that 
for "typical" realizations (x\,- • ;x m ), we have 

n 

— = P 
m 

and, hence, the two distributions have the same general location and 
scale. [In fact, expression (31) is a hypergeometric distribution, which 
converges to (28) as m — - °° if n = pm (Ref. 13).] Thus, although incor- 
rect, the binomial distribution approximates the true distribution of the 
killer trunk. 

The following result helps to put the relationship between the ad hoc 
and the optimal statistics in perspective. 

Lemma 2: If |x,} is a binary state stationary Markov chain with tran- 
sition probabilities P ,i and Pi , and if x m = {xi,--;x m ), then we 



1 Pi.oip.r) —■ 1 - p as r — °°. 
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have 

P(x m ) = b(ti ;n,Pi,o) X b(t ,i;m-n,P 0il ) X q, (33a) 

where 



P(si)(l-Pi.o)-*"(l-Po,i)- 

q = 



(33b) 



C )( m ~ n ) 

\tio/\ t i I 
Proof: 

P( Xm )=P( Xl ) UP(Xi/Xi-l) 

i=2 

— P(y,\ P' 10 P'" P' 01 P' 00 

— r\x\) r lj0 r xx r 0A r Q0 

= PQ (1 - P,.o)'" P l a ! (1 - Po.i)' 00 X P(xi). 

Using eqs. (11a) and (18) to express tn and too in terms of iio and £ i, 
respectively, yields the result. 

Thus, given a binary state stationary Markov chain [*,•}, it is clear from 
the above lemma that the log likelihood ratio 

i(x m ) = log-^ - 

formulated for the two hypotheses 

H : \xi\ Markovian, characterized by (Po.i.Pi.o) 
ffi: \xj\ Markovian, characterized by (Po,i,Pi,o) 
is the sum of three terms: 

*/ v , fe(tio;n.pj.o) . , b(t i;m - n,Pl.i) . /q*\ 

e(x m ) = log— ^-r+ log 77: E^ + log (T/- 

b(t l0 ;n,P h0 ) b(toi;m - n,P ,i) \ q I 

The first term is the ad hoc statistic (at 10 - an), the second term is the 
additional statistic [/ftoi — b(m — n)], and the third term is an end-effect 
term (ax m + bx l m ). 



log (—) =ax m +bx l m + log 



P*(xi) 
P(xi) 



= ax m + bx 



rm 



since,P*U,) = PUi) = p. 

The ad hoc algorithm, although based on the approximate binomial 
distribution, is very attractive for a number of practical reasons: 

(£) For the 200-second sampling option, it is essentially optimum in 
a practical sense, since the P0.1 characteristics for a normal and killer 
trunk are not far enough apart to exploit (see Fig. 6). 

(ii) When we exploit the grouping information for the 5XB trunk 
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group in Section V, it will become obvious that the additional part of the 
optimum statistic is quite sensitive to the trunk occupancy and, hence, 
to the trunk-selection procedure modeled. (The sensitivity of the "ad- 
ditional" statistic to trunk occupancy stems from the fact that Po.i is 
"almost" proportional to p.) On the other hand, we will see that the ad 
hoc 5XB algorithm is relatively robust to minor perturbations in the 
trunk occupancy (and therefore to the trunk-selection procedure) and 
hence might be expected to perform well in a real 5XB environment. 

V. THE 5XB TRUNK-GROUP ALGORITHMS 

In addition to utilizing individual trunk switch-count and transition 
accumulations the 5XB group algorithms exploit the following: 

(0 The identity of all trunks common to a group. 
(ii) The trunk-selection procedure. 

The resulting 5XB group algorithms typically are faster* than their 
individual trunk counterparts and are also less sensitive to the groups 
nominal holding time. 

5. 1 The 5XB trunk-group model 

For the purposes of this paper, we model a 5XB trunk group (with all 
trunks normal) as an M/M/N-loss model with random selection of idle 
trunks. 2 The same assumptions apply if the group contains one or more 
killer trunks, but in this case we assume that killer trunks have a mean 
holding time equal to l/r that of the normal mean holding time. In ad- 
dition to being convenient theoretically, this idealized model has also 
been very useful in developing the 5XB group algorithms presently 
implemented in ICAN. 

If all N trunks are normal, the random selection rule implies that all 
trunks have the same mean occupancy. In Ref. 2, the birth and death 
equations for the above model with a single killer trunk were solved in 
closed form, and in Ref. 14 this was generalized to an arbitrary number 
of killer trunks. These analytic results turn out to be quite useful, and 
in what follows we will need the following results derived in Ref. 2. 

Theorem 1: For the above 5XB trunk group model having a single killer 
trunk with parameter r and an offered load of a erlangs, the blocking 
probability &{N,a,r) and the mean occupancy p*(N,a,r) of the killer 
trunk are given by: 



* For given type 1 and 2 errors, the group algorithms typically have a considerably 
smaller mean decision time than their individual trunk algorithm counterparts. 
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(») P ;(N,a,r) = , (34b) 

1 + — -r[l-B(N-l,a)] 

a 

where B(N,a) is the usual erlang B blocking associated with an Ml 
M/N-loss system with all trunks normal and an offered load of a 
erlangs. 

It is easy to see that the occupancy p r of each of the N - \ normal 
trunks must satisfy the conservation equation: 

rp* r +(N-l) Pr = a[l-&], (34c) 

and the trunk-group occupancy 0* is defined by: 

,, P ; + (N-i)p r 



N 



(34d) 



(For a 5XB trunk group having a killer trunk with parameter r: p* and 
p r denote the mean occupancy for a killer and normal trunk and 0* de- 
notes the mean group occupancy.) 

Although eqs. (34a) through (34d) define an implicit relationship 
between p*{N,a,r) and 0*(N,a,r), it will be very useful to have a simple 
explicit relationship. If the blocking term in eq. (34b) is ignored and if 
we "associate" 0* with a/N, an approximation suggested is: 

(35a) 



"' r-(r-l)0r 

This approximation, although quite good for large N, is rendered ob- 
solete by the following exact result: 

Theorem 2: Consider a 5XB trunk-group model with all trunks normal 
and mean-group occupancy 0. (We will let denote the (mean) group 
occupancy for a 5XB trunk group with all trunks normal.) If one of the 
trunks is replaced by a killer with parameter r, then 

P* = p(0,r), 



where, 



Of course 



P(0,r)= 7 -- . (35b) 

r - (r - 1)0 



q[l- B(N,a)] 
N 

TRUNK-DETECTION ALGORITHMS 947 



where B(- , •) is the usual erlang B blocking expression. 

This surprising result, which follows easily from eq. (34b), is proved 
in Appendix D. As a consequence of this theorem, "5XB group occu- 
pancy" will be used to denote the occupancy of a 5XB group with all 
trunks normal. 

The mean occupancy of the normal trunks in a 5XB trunk group 
model having a single killer trunk no longer is given exactly by 0. But 
the following result, derived in Appendix D, shows that is a good ap- 
proximation. 

Theorem 3: Consider a 5XB trunk group model with N trunks having 
a single killer trunk with parameter r ^ 1. The mean occupancy (p r ) 
of the N — 1 normal trunks satisfy 

N 
r 



-(*£>-»*; 



*Pr*l : — JX0> (36) 

r - (r - l)4> 



where <f> is the mean-group occupancy with all trunks normal. 

Theorems 2 and 3 are proved in Appendix D, where an exact expression 
for p r is also derived. These results are special cases of general results 
obtained for the random selection model. 14 

5.2 Exploiting the 5XB Grouping Information 

To simplify matters, we assume that a trunk in a 5XB group* with 
mean-group occupancy <f> has mean occupancy p(<f),r) given by 

P(<t>,r) = / — , (37a) 

r - (r - 1)0 

where r = 1 corresponds to a normal trunk. Thus, if a group has no killer 
trunks, all normal trunks satisfy p = <f> and eq. (37a) with r = 1 yields the 
correct occupancy. If, however, the group has a killer trunk, then all 
normal trunks satisfy inequality (36) so eq. (37a) with r = 1 is an ap- 
proximation that increases in accuracy with the size of the group. Of 
course eq. (37a) is exact for a (single) killer trunk in a 5XB trunk 
group. 

It is clear from eq. (37a) that 

"^ m f P< *nui> ,37b) 

r - {r - l)p(0,l) 



* We use "5XB group" and our idealized model of a 5XB trunk group interchange- 
ably. 
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and hence p{<b,r) is typically much smaller than p(#,l). (For <f> = 0.50 and 
r = 10, p((j),r) = 2/11 p(0,l). Thus, it would appear that considering the 
1 — *■ transition probability as a function of the 5XB group occupancy 
would effectively "spread" the Pi,o characteristics in Figs. 2 and 3 further 
apart. That is, for a given 0, we propose comparing 

Pi,o[p(<M)] and P liO [p(0,r)] 

[rather than Pi,o(p) and Pi,o(p), as in Section IV]. 

Denoting the composition Pi,o[p(0,^)] by P\,o{<t>,r), we have 



-rS 



P 1 ,o(0,r) = (l-p(0,r))(l-expj T3 ^)), 



(38) 



which is plotted in Figs. 8 and 9 for the 200- and 100-second sampling 
options, respectively. The normal holding time used in these figures is 
180 seconds, and the killer-trunk characteristics are drawn for r = 5, 10, 
and 15. 

The increased "spread" between normal and killer Pio characteristics 
obtained in this way is simply a consequence of exploiting the distinctly 
different occupancies of a normal and killer trunk in a 5XB trunk group. 
Figure 10 is a three-dimensional sketch of the composition of Pi,o and 
p. Because all normal trunks in a 5XB group have the same mean occu- 
pancy, we see that a single Pi vs 4> characteristic suffices to describe 




5XB TRUNK GROUP MODEL 
(RANDOM SELECTION OF IDLE TRUNKS) 
-' = 180 SECONDS 
T= 200 SECONDS 



20 40 60 80 

MEAN GROUP OCCUPANCY IN PERCENT 



100 



Fig. 8 — 1 -* transition probability for the 200-second sampling option. 
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5XB TRUNK GROUP MODEL 
(RANDOM SELECTION OF IDLE TRUNKS) 
/i"' = 180 SECONDS 
T= 100 SECONDS 



20 40 60 80 

MEAN GROUP OCCUPANCY IN PERCENT 



100 



Fig. 9 — 1 -* transition probability for the 100-second sampling option. 




P 10 .pPLANE:-- 

M/M/1 - LOSS 
MODEL 



7 Z p.0PLANE: 5XB GROUP 

/ MODEL 



p(0.r 



Fig. 10 — Sketch of the composition of Pi,o(/>) with p(<£,r). 

all normal trunks. This fact allows us to translate the individual trunk 
algorithm's development to this 5XB context with essentially only no- 
tational changes. 
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20 40 60 80 

MEAN GROUPOCCUPANCY IN PERCENT 



100 



Fig. 11 — 1 -► transition probability for the200-second sampling option (mean group 
holding time = 45 seconds). 

Figures 11 and 12 are plots of eq. (38) drawn for a normal group mean 
holding time of 45 seconds.* The normal trunk characteristic corre- 
sponding to 180 seconds is shown in dashed lines. We see that with 5XB 
grouping information factored into the picture, considerable discrimi- 
nation exists between both normal trunk characteristics as well as be- 
tween the normal trunk having a holding time of 45 seconds and the killer 
trunks. The discrimination that exists between the normal trunks per- 
mits us to make the 5XB group algorithm adaptive to the group mean 
holding time. [Although we will not pursue this topic, the basic idea is 
that Ey£io(/)/2/n(/) (sums are over all trunks in the group) is an estimate 
of p i,o and can be used to decide which (of several) normal pi,o charac- 
teristics constitutes Ho.] 

5.3 The ad hoc and the optimal 5XB group algorithms 

We assume that the mean group holding time is known and consider 
formulating a hypothesis-testing problem similar to that in Section 4.1. 
Thus, we denote P(t iQ (m) = t/n(m) = n) by P{t/n) and consider the two 
hypotheses: 



* For normal holding times in the vicinity of 45 seconds, a killer parameter r in the range 
3 to 5 probably is typical. An r of 10 or 15 in this context is unrealistic. 
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20 40 60 80 

MEAN GROUP OCCUPANCY IN PERCENT 



100 



Fig. 12 — 1 — transition probability for the 100-second sampling option (mean grouping 
holding time = 45 seconds). 



H : P(t/n) = b[t;n,P h0 (<t>,l)] 

Hi: P(t/n) = b[t;n,P h0 (<i>,r)], r £ rn. 



(39) 



There are two differences between this formulation and the one in 
Section 4.1: 

(i) The trunk occupancy p in Section 4.1 is replaced by the group 
occupancy <f>. 

(ii) The alternate hypothesis H\ is composite since P\,o(<t>,r) for r ^ 
ro are distinct. 

The approach taken in dealing with (ii) is a natural one often 
adopted; 11 since Pi, o(<t>,r) is monotone increasing in r {this follows from 
eq. (38) upon noting that r/[l — p(0,r)] = r + 0/(1 — <J>)\, then testing 
between Ho and the simple alternate hypothesis 

HnP(t/n) = b\t;n,P lt o{4>,ro)l 

say with type 1 and 2 errors a and /3, respectively, implies that if the true 
state of nature is Hi with r = r > ro the resulting type 2 error will not 
exceed 0. With this approach, we can simply translate the ad hoc algo- 
rithm results developed in Section 4.1 to this 5XB group context by 
making the appropriate changes in notation. 
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Thus, the ad hoc 5XB group algorithm can be described as follows: 
corresponding to the sequence of data accumulations (t\,n\ ) i - 1,2,- • • 
for the jf th trunk in a group of N trunks, where t{ and n{ are the 1 -* 
transition and switch-count accumulations, respectively, during the ith 
accumulation period in which m, scans are made, define s\ and S\ i = 
l,2,...;; = l,...;Nby 

si = a(4>i)ti - a(<t>i)nj (40a) 

and 

Si = Si- X + si with Si = 0, (40b) 

where <fr is the group's occupancy during the ith accumulation period. 
The sequential test for the ;'th trunk in the group, ;' = 1,- • -JN is defined 
by 

(0 Compute S{, i = 1,2,* • • and defer making a decision as long as T 
< Si <T h 

(ii) Ifi = k corresponds to the first accumulation period for which 

Si * {T Q ,T X ), 
then 

S{ ^ To =* trunk ; is normal 

S{ ^ T\ =* trunk ; is a killer. 
The weights a(0) and «(</>) are defined by 

fl( * , - l0g l-Pt.aC0.ro) (4U) 

and 

«(*)-«W + log£^, (41b) 

where Pi,o(<l>,r) is defined by eq. (38). 

Just as in the individual trunk algorithm, the actual occupancy re- 
quired to choose the weights a and a is unknown and must be estimated. 
Thus, the group occupancy tf>, during the ith accumulation period is 
estimated by 

ii = h £ "< /m < (42) 

Note that fa is a "better" estimator than p, = n,/m, (the estimator used 
in the individual trunk algorithm) in the following sense: given a 5XB 
group with all trunks normal and mean-group occupancy (in equilib- 
rium), we have 
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(l) E(p) = p = 

(w) E(0) = 

(m) var(0) < var(p). 

In addition to the better occupancy estimate available on a group basis, 
the fact that the group Pi$ characteristics are "flatter and broader" then 
the individual trunk Pio characteristic implies that the group algorithm 
more faithfully tracks the required weights, than does the individual 
trunk algorithm. 

The ad hoc 5XB group algorithm has the same pleasant intuitive in- 
terpretation that the ad hoc individual trunk algorithm had (see eq. (13) 
and related discussion). It is also easy to show how the optimal individual 
trunk algorithm development of Section 4.2 carries over to the 5XB 
group context. 

Thus, consider the two states of a trunk to be described by: 

Ho: \xi\ Markovian, characterized by = (Po,i,Pi,o) 
H\. \xi\ Markovian, characterized by 6* = (Po,i,Pi,o)» 
where Pi.o = Pi.oOM), Pi.o = Pi.oi^o), and [see eq. (2)] 

PoMr) = - P(0 ; r j , Pi.o(*,r) (43) 

1 - p(0,r) 

with p((/>,r) defined by eq. (37a). The assumptions that lead to a con- 
sideration of these two statistical hypotheses as a model of the normal 
and killer states of a trunk can be found in Section 4.2. 

Proceeding as in Section 4.2 leads us to the optimum statistic £(x m ) 
for distinguishing between the two simple hypotheses under consider- 
ation: 

Hx m ) = [(a - a)t 10 - at n ] + [Bt 00 -0- B)t 01 ] + log ^7^ , (44) 

PUi) 

where the parameters a and a are defined by eqs. (41a) and (41b) and 
the parameters b and are defined by 

r , 1 ~ Pq, i(0,ro) ... x 

6 = log i-p„,<*,i> (45a) 

and 

P = b + \og • (45b) 

Po,i(4>,ro) 

As in eq. (39), the alternate hypothesis is really composite since 
[Po.ii^yf), Pi,o(<t>,r)] for r ^ ro are distinct. Since Po,i(4>,r) is monotone 
decreasing in r, the approach discussed earlier of treating H\ as a simple 
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hypothesis with r = r is followed. The parameters b and p are defined 
a bit differently than were the parameters b and for the optimal indi- 
vidual trunk algorithm (see eq. (20a) and (20b)) in order to obtain non- 
negative weights. Thus, in the individual trunk algorithm context, we 
had Po.i > Po.i (eq. (19c)), but in the present 5XB group context we have 
Po.i > P'o.i- Tne reason for this "flip-flop" is easy to see: for a given group 
occupancy <t>, we are now contrasting the 0—1 transition probability 
for two trunks which differ, not only in their hang-up rates but also in 
their occupancies as well. Thus, the difference in the two occupancies 
dominates the effect that the hang-up rates alone have. Roughly 
speaking, the — 1 transition probability of a trunk is approximately 
equal to its occupancy (conditioning on the last scan has little effect) 
and hence since p(^,r Q ) « p(<£,l) it is clear that we should have Po.i(<£,r) 
< P ,i((f>,l)- Figures 13 and 14 are plots of Po,i(<t>,r) for the 200- and 
100-second sampling options, respectively. In both figures, the killer 
trunk characteristics have been plotted for r = 5, 10, and 15 and a normal 
mean holding time of 180 seconds is assumed. These figures are very 
insensitive to the assumed normal mean holding time, since they es- 
sentially reflect eq. (37a), which is independent of the mean group 
holding time. 

The "additional" statistic which appears in eq. (44), 

Bt oo - - 6)t i 
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Fig. 13— -* 1 transition probability for the 200-second sampling option. 
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MEAN GROUP OCCUPANCY IN PERCENT 
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Fig. 14 — -* 1 transition probability for the 100-second sampling option. 



shows that — * transitions are weighted positively (evidence of a killer 
trunk*) and -* 1 transitions are weighted negatively (evidence of a 
normal trunk). This additional statistic is strongly influenced by the 
occupancy of a trunk, and only slightly by its hang-up rate. 

Note also that the term log P* (xi)/P(x i) is nonzero in the 5XB context 
since 



if xj = 1 
if xi = 0. 





1 


P*(xi) 
P(xi) 


\r-(r-l)<t> 
) r 




V-(r-l)0 



VI. PERFORMANCE OF THE 5XB GROUP ALGORITHMS 

In common with all sequential detection algorithms, the time required 
by the killer-trunk detection algorithms to reach a decision (trunk nor- 
mal or killer) is a random variable. In this section, we obtain an ap- 
proximate formula for the mean time required by the 5XB group algo- 
rithms to reach a decision. This result is used to contrast the performance 



* Killer trunks in the 5XB group model have very low occupancy, and hence — * 
transitions are likely. 
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of the ad hoc and optimal algorithms as well as to point out the consid- 
erable effect that the sampling rate has on each. In addition, we also 
obtain an approximate expression for the false-alarm probability of the 
5XB group algorithms. The analysis for the individual trunk algorithms, 
although differing in several respects from the group algorithms, involves 
the same sort of considerations and is omitted. 

The analysis in this section assumes a server system in equilibrium 
and, therefore, the mean trunk-group occupancy is assumed constant. 
In addition, to simplify the analysis, we assume that is known; an ap- 
proximate analysis which does not require this assumption is sketched 
in Section 6.1. A consequence of this assumption is that the algorithm 
weights are treated as constants rather than random variables. This 
assumption is not unreasonable because for multiple-hour accumulation 
periods var(0) is quite small (0 is the switch-count estimate of 0). (Var(0) 
has been derived for an M/M/N-loss system. 15 ) 

6. 1 Mean statistic update 

Corresponding to a sequence of trunk states x x ,x 2 ,- • -,x m in an accu- 
mulation period with m scans, define a sequence of transition updates 
z>i,- • ;z m by 

6 if (x»-b*») = (0,0) 
( - (0-6) if (x„_i,x n ) - (0,1) (46) 

|(a-a) if (x„-i,x„) - (1,0) 
-a if (Xn-ifXn) = (1>1) 
The optimum 5XB statistic (eq. (44)) may therefore be written: 

m P*(x m ) 

/(x m )=L^ + log^T (47) 

, = 2 P(x m ) 

In practice the end-effect term cannot be implemented and all the 
transitions in eq. (44) must be estimated in terms of t{m), n(m), and m. 
Thus, if we denote the implementable version of eq. (44) by S m (0,r o ), t 
use eqs. (11), (12), and (20) in eq. (44), and drop all end-effect terms we 
obtain 

S m (0,r o ) = «(0) — a((p)n(m) 

+ |^(0)[m-n(m)]- 0(0)^1 (48) 

f r is the value of the killer parameter used in defining the alternate hypothesis and 
hence the algorithm weights (see eqs. (41) and (45)). 
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Zn = 



Now using eqs. (11), (12), and (20) once again, we see that S m (0,ro) may 
be written 



S m (<f>,r ) = £ Zi + (bx c m - ax m ) + (« + /&) l^— ^). 

1-2 \ Z / 



(49) 



Since the 5XB group is in equilibrium, z%- • -,z m are identically dis- 
tributed and it is easily verified that their common mean is given by 



E(z) = (a Pi, - a)p + (b - P o ,i0)(l - p). 



(50) 



We recall that p (eq. (37a)) as well as the transition probabilities (eqs. 
(38) and (43)) are functions of the group occupancy and the state r of 
the trunk. The weights a, a, b and j3 are only functions of for a specified 
choice of the parameter (ro), which characterizes the alternate hypoth- 
esis. Thus, the mean statistic update for the optimum 5XB statistic and 
its "implementable version" is given by: 



E\kx m )\= (m-l)E(z)+E 



log 



P*(*i) 



P(*i) 



and 



E[S m (<p,r )\ = (m - l)E(z) + 5(1 - p) - ap. 
Note that it is easily shown that 



(51a) 



(51b) 



E 



lo g r,, v ■ -log [ro - (r - 1)0] + 



P(x : 



r-(r- 1)0 



logr 



which is negative for r = 1 and positive for r ^ ro- 

Although the increments (transition updates) defined in eq. (46) are 
identically distributed, they are not independent. In fact, since the state 
sequence \x, \ has been modeled as a Markov chain (see Section 4.2), it 
is easy to see that the sequence {z,-} defined by (46) is also a Markov chain. 
Relabeling the four natural states, 



6, -(0-fi), (a -a), -a, 



of the chain {z,-} by 









0, 1, 2, 


3, 








respectively, 


it is easily 


seen that th€ 


one- 


step transition matrix ir for 


this chain is 


given by 
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(a) 

OPTIMAL STATISTIC 



1-Pi n> 




(b) 

AD HOC STATISTIC 

Fig. 15 — State diagram for the ad hoc and optimal 5XB group statistics, (a) Optimal 
statistic, (b) Ad hoc statistic. 



where ir,j = P(z n = }lz n -\ = i). Of course, the stationary distribution 
P satisfying Pit = P is 

P = [(1 - p)(l - Po.i), (1 " p)Po.l pPi.o, p(1 - Pi.o)]- 

The state diagram corresponding to the Markov chain (z, J is shown in 
Fig. 15a. If b = = 0, we obtain the ad hoc algorithm, for which the se- 
quence \zi\ is a three-state Markov chain, with natural states 0, a — a 
and —a. The state diagram for this chain is shown in Fig. 15b. 

Figures 16 and 17 are plots of the mean statistic update vs group oc- 
cupancy for the implementable version of the 5XB group algorithms. 
(Log base 10 is used in this paper.) These figures are drawn for a killer 
trunk with parameter r = 10. Also shown is the corresponding plot for 
a normal trunk (r = 1). 

We close this section by indicating an approximate analysis of E{z) 
which doesn't assume that <t> is known. Thus, in practice, is unknown 
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Fig. 16 — Mean statistic update for the 5XB group algorithms — 200-second sampling 
option. 
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Fig. 17— Mean statistic update for the 5XB group algorithms — 100-second sampling 
option. 
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and is estimated by 0. Hence eq. (46) should read: 

(6(0o) if (x n -i,x n ) = (0,0) and = O . 

Conditioning on 0, we have 

E(z) = E\E(z/b)\, 

where 

E(z/0) = 6(0)P(x„-i = 0,x„ = 0/0) + • • .. 

Now by assuming that (x n -i,x n ) and are independent,* we get eq. 
(50) with the constant weights «(0),- • • replaced by the mean values 
E|a(0)},- • •. The mean values can be approximated in either one of two 
ways: 

(i) £ja(0)j = «[£(0)] = rt(0), which obviously amounts to assuming 
is known. 

var(0) d 2 « t 



(ii) E\a(<f>)\ = «(0) + 



2! d0 2 
which factors the available variance of into the picture. 

6.2 Mean time to detection 

The basic structure of all the detection algorithms in this paper are 
the same: a statistic s, is evaluated at the end of the ith accumulation 
period, i = 1,2,- • • and a decision is made the first time that the sum Si 

+ S2 H falls outside an interval (To,Ti). Presumably, the random walk 

type statistic S, has a negative drift under Ho (trunk normal) and a 
positive drift under H\ (trunk killer). Wald's SPRT always has the ap- 
propriate drift: if Ho and H\ correspond to the probability distributions 
Pq(w) and P\(w), respectively, and if Hw) is defined by 

Po(w) 

then E\l{w)\ < under H and E{i(w)\ > under H\. The proof is im- 
mediate by using the inequality 12 

- HP, log Pi < - £ Pi log q i} 



* For reasonable-size trunk groups, we expect very little dependence between the 
sampled state process of an individual trunk (.r i,- • •) and the group process <Jt. 
f This is a Taylor expansion to second order. 
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where {p,| and \q, ] are distinct probability distributions. In Appendix 
E, we show that both the ad hoc (6 = ft = 0) and the additional (a = a 
= 0) parts of E(z) (eq. 50) have the appropriate drift. This in turn shows 
that both the ad hoc and additional parts of the 5XB group statistic (eq. 
(51a)) have the appropriate drift. 

Suppose Y\,Y2,' ' ' are i.i.d. random variables with common mean n 
and consider a random walk 

Sn=£.Y i n= 1,2,- • • 

with absorbing barriers at To and T\. If n is small compared to To and 
Ti, then the mean stopping time (mean number of steps to absorption) 
E(n) is approximately given by 

PoTo + PiTi 



where Pq and Pi are the probabilities of absorption at To and T\, re- 
spectively. This follows from Wald's identity 11 

E(S n ) = nE(n) 

if we approximate the mean value of the random walk at absorption by 
P T + PiTl 

In our detection theory context, Po and Pi correspond to /3 and 1 — 
/?, respectively, if the trunk is a killer (/? = probability of miss) and 1 — 
a and a, respectively, if the trunk is normal (a = probability of false 
alarm). If we denote the mean number of accumulation periods needed 
to reach a decision under Ho and H\ by E(Tn) and E(Tk), respectively, 
and assume (i) successive statistic updates are independent and (ii) the 
mean statistic update is small compared to T and T\, we obtain 

E{Tn) = CiSUtetH (52a) 

and 

. /3T +(l-/3)Ti 

EiTk) - E\S m (*,r )\ ' (52b) 

where the mean statistic update E\S m (0,ro)j is evaluated for r = 1 in (52a) 
and for r s ro in (52b). 

The assumption that successive statistic updates are independent is 
not strictly true if the successive statistic updates are contiguous. 
However, one expects that the slight (end-effect) dependence will not 
give rise to very much error. 

Let Tk be the time required to decide (incorrectly) that a killer trunk 
is normal. Similarly, let T^ , be the time required to decide (correctly) 
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that a killer trunk is a killer. Then we have 

E(T k )>*0E(T ko ) + {l-(3)E(T kl ) f 

which suggests that 

E(T ki )±E(T k ) if 0«1. (53) 

The moments of the conditional stopping times Tk and Tk , can be ob- 
tained by using a well known technique of Wald's,* and one finds that 
approximation (53) is reasonable if T\ and (—To) are sufficiently 
large. 

In our faulty-trunk detection context, a type 1 error ("false alarm") 
may result in the misuse of craft resources (e.g., testing a perfectly good 
trunk). A type 2 error ("miss") on the other hand, will result in an in- 
creased time to detection. Assuming type 1 and 2 errors of 10 -6 and 10~ 2 , 
respectively (realistic implementation values), implies approximate 
thresholds T = —2 and T\ = 6 (formulae 4.5c and d, log base 10). With 
these parameter values, the mean detection time E(Tk x ) can be ap- 
proximated by 

E(T "^m^y (54) 

For a normal trunk, the mean statistic update is comparable to To and 
hence expression (52a) isn't applicable, nor is it needed since the 5XB 
group algorithm reaches a decision on a normal trunk after one or two 
updates. 

Figures 18 and 19 are plots of the mean-detection time vs group oc- 
cupancy for the implementable versions of the 5XB group algorithms. 
[The dashed line portion of Figs. 18 and 19 indicates where approxi- 
mation (54) involves considerable error (e.g., the region in which 
E\S m {<j),r{))\ is a significant fraction of T\).] It is apparent from these 
figures that 

(0 The mean detection time for both the ad hoc and optimal algo- 
rithms is enhanced by using the 100-second sampling option. This en- 
hancement is far more pronounced for the ad hoc algorithm. 

(ii) The optimal algorithm is "faster" than the ad hoc algorithm. This 
contrast is greater for the 200-second sampling option. 

6.3 False alarm probability of the 5XB group algorithms 

If B{ot,fl) and A (a, (3) are the test thresholds that result in type 1 and 
2 errors a and 0, respectively, for a SPRT, Wald showed 11 that using 



* See eq. 158 and 159 in Appendix A. 5.2 of Ref. 11. The method ignores (as usual) the 
'excess over the houndaries" and hence yields approximate results. 
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Fig. 18— Mean time to detection for the 5XB group algorithms — 200-second sampling 
option. 



thresholds To and T\ defined by* 



To = log (tz~) and Tl = log (~/ ( assumin & « + < D 



yield type 1 and 2 errors a' and 0' which satisfy 



0'* 



and 



l-« "1-/? 

The proof of this result is trivial and depends only on the assumption 
that the SPRT terminates with probability 1. This assumption is satisfied 
by a wide class of SPRTs, 16 including the case of interest to us, where the 
underlying distribution is "Markovian". Thus, the probability of false 
alarm («') for a SPRT (an example of which is the "optimal" 5XB group 
algorithm) satisfies «' ^ 10" T K 

Because the ad hoc algorithm is not a SPRT (the true underlying dis- 
tribution is not binomial, see Section 4.3), we may wish to study the 
consequences of using the above thresholds appropriate for a SPRT. To 
do this, we can proceed in (at least) two distinct ways: 

(i) consider the statistic S m H,ro) to be the basic update in the algo- 
rithm. 



* This assumes that the log likelihood ratio is used in defining the SPRT. 
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Fig. 19 — Mean time to detection for the 5XB group algorithms — 100-second sampling 
option. 



(ii) Consider the transition update Z, (x,_i,x,) to be the basic update 
in the algorithm. 

In either case, we study a process of the form Sj = 2f =1 xi, j = 1,2, — 
with absorbing barriers To < < T\. The advantage of proceeding as in 
(i) above is that the increments x, may be assumed to be i.i.d. We briefly 
sketch this approach. 

Consider a random walk Sj, j = 1,2,« • • with i.i.d. increments {#;}. If 
n is a stopping time associated with Sj, j = 1,2,- • •, Wald's fundamental 
identity 1117 is given by 



where 



E\e s "'x(t)~ n \ = 1 (all t satisfying \xit)\ ^ 1), 



x(t) = E\e*t] 



(55) 



is the moment generating function corresponding to the common dis- 
tribution of the increments. If Pt q and Py, denote the probabilities of 
absorption at To and T\, respectively, then rewriting (55) in a standard 
way yields 

P T{) E\e»»<x(t)-»/S n 5; T | + P Tl E\e s '> t x(t)- n /S n * T,| = 1. (56) 

Now if the x, take on both positive as well as negative values with nonzero 
probability and have a nonzero mean, then the equation 

x(0 = i 
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has only one nonzero root /in with the property that ho and E(x) have 
opposite sign.* That is, E(x) > =» h < and E(x) < =» h > 0. As- 
suming that the excess of S„ over the boundaries is small, eq. (56) yields 
the standard approximation 1117 

p hoTi _ i 

P To * Ti; I^r (57a) 

e hoTi _ e hoTo 



and 



i _ p/ioTo 
1 e noTi _ e hoTo 



Note that the probability of false alarm (type 1 error) corresponds to Pti 
if the random walk increment is that of a normal trunk [S m (0,ro), r = 
1]. Similarly, the probability of miss (type 2 error) corresponds to Pt 
if the random walk increment is that of a killer trunk [S m (0,ro), r ^ 

roL 

To use these approximations, we must compute the moment gener- 
ating function 

X (») = E\e M*)«/2-a (0)n |m} (58) 

by using the joint distribution p{t,n) derived in Appendix B. [Note that 
our discussion applies equally well to the additional and optimal sta- 
tistics. In general, we need E\exp [S m (0,r o )ju]l, where S m (</>,r ) is given 
by eq. (48).] Choosing the test thresholds To and T\ for the ad hoc al- 
gorithm according to the Wald SPRT formulae [eqs. (10c) and (10d)] 
typically results in Pr < fi and Pt x < «• 

VII. SUMMARY 

A class of killer-trunk detection algorithms has been developed that 
use the individual trunk usage and transition accumulations available 
in EADAS/ICUR. Because this data is essentially a sufficient statistic for 
the Markov chain used to model the (unobservable) sampled data, one 
of the algorithms developed is Wald's celebrated SPRT. 

The detection algorithms developed can be partitioned in two natural 
ways: 

(i) By sampling rate (100 or 200 seconds). 

(ii) According to whether grouping information is used. 

The algorithms which do not use grouping information are applicable 
to all trunks (one way or two way) independent of the type of switching 
machine used. A version of one of these individual trunk algorithms is 
currently in use in ICAN, testing trunks on 1XB, XBT and step-by-step 



* See Appendix A.2.1 of Wald's original treatise (Ref. 11). 
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switching machines. The algorithms that exploit grouping information 
detect killers more quickly but are tailored to a specific switching ma- 
chine. A "group" algorithm of this type is currently being used to test 
trunks associated with 5XB switching machines. 

In the course of this study several problems of independent interest 
were studied. These include: 

(i) The server covariance in a M/G/1-loss and GI/M/1-loss system. 

(ii) The structure of the likelihood statistic that arises in testing 
simple hypotheses characterized by a binary valued Markov chain. 

(Hi) The occupancy of nonidentical trunks in a random-selection 
(Markovian) loss system. 

The major conclusion in this study is, of course, that accumulated 
switch-count and state-transition data on individual trunks (based on 
sampling intervals on the order of a normal holding time) can be used 
to reliably detect abnormally short holding time trunks. Moreover the 
(near) optimal sequential detection algorithms using this accumulated 
data are easily exhibited, simple in structure, and intuitively appeal- 
ing. 
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APPENDIX A 

Sensitivity of the Transition Probabilities to Modeling Assumptions 

To get an idea about the sensitivity of the algorithms to some of the 
modeling assumptions, the transition probability Pi r o = P(x t + T = 0/xt 
= 1) was studied for the following two cases: 

(i) M/G/1-loss, where the service distribution function F(-) is the 
mixed exponential given by 

F(t) = 1 -die-*"' -d 2 e _X2 ' i^O 

(ii) GI/M/1-loss, where the arrival process is the switched-Poisson 
process 8 commonly used to model overflow traffic. 

Because the methods used to obtain Pio for these two models differ, 
we discuss these models separately. 

A.1 The M/G/1 -loss model 

An observer viewing the server in an M/G/1-loss system sees an al- 
ternating sequence of busy and idle intervals. The busy intervals are 
distributed according to some distribution F(-) and are independent. 
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x t = 



The idle intervals are exponentially distributed with mean X -1 (X = mean 
arrival rate) and are independent. Thus, the sequence of alternating busy 
and idle intervals constitutes an alternating renewal process. In this 
context, the conditional probability 

Pi.o = P(x t +r = 0/x t = 1), (59) 

where 

1 if server is busy at epoch t 

if server is idle at epoch t 
has already been studied, 18 and we have the following result: 

Theorem 4: Consider an M/G/1-loss system in equilibrium with a ser- 
vice time distribution F(t) having Laplace transform f*(s). If P\ (s) 
denotes the Laplace transform of P\q{t), then we have 

s\s+ X[l -f*(s)\\ 

where X and n are the mean arrival and service rates, respectively. 

Proof: See Section 7.4 of Ref. 18 [f* 2 (s) = X/s + X and f\(s) = /*(«)]. 

For the mixed exponential service time distribution mentioned above 
let 

xi = 2nd 
and 

x 2 = 2n(l-d), O^d^l. 
If T is distributed according to this two-parameter (n,d) family of dis- 
tributions, then 

(0 E{T) = >l-K 

(ii) var(T) = m -2 (^T - ). where « - 4d(l - d). 

,-.-x ,™ °(T) /2 - 5\ 1/2 

u " um =wr(— ) • 

Thus, the mean is fixed at n~ l and the coefficient of variation satisfies 
c 3: 1, with equality occurring when 5 = 1, which corresponds to the 
M/M/1-loss system. 

Equation (60) is easily inverted for this family of mixed exponential 
distributions, obtaining the following result: 

Pi.o(p.r) = (1 - P )+ n + 5 - g nMT+ r2 + 5 e ^r (61) 
r\(ri-r 2 ) r 2 {r2-r 1 ) 

where 
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r\ = 



r-> = - 



f-±} + [a -«<i +.> + ©■] 

)-[ (1 -a> (1 + fl) + ( |)°] 



1/2 



2 

b = Ad(l-d), a=—^~ 
1 - p 

and p = X/(A + p) is the mean occupancy of the server. 

In Fig. 5, P\,{)(p,t) [eq. (61)] is plotted vs p for several different values 
of c (c = 1, 1.5, and 2) assuming 

(i) p~ l = 180 seconds 
{ii) t = 200 seconds. 

Also shown is a plot of Pi, () vsp for a killer trunk with r = 10 [p. in (61) 
is replaced by rp]. The normal trunk Pi, characteristic with c = 1 cor- 
responds to the M/M/1-loss system. 

A.2 The GI/M/1-loss model 

The covariance function /?(•) for the GI/M/1-loss model with a 
switched Poisson arrival process has the form: 

R(t) = c, e w ' r + c,e w * T + c 3 e- {w+ y )T , (62) 

where the coefficients c, and the exponents Wi are messy expressions 
involving the three switch parameters co, y, and X (Ref. 8) and the mean 
service rate p. The derivation of this covariance function is straightfor- 
ward but tedious and is therefore omitted. (For the switched Poisson 
arrival prccess, the Markovian state equations can be solved for 
P(x,,Xi+ T ), where x t = state of server at epoch t.) Our purpose here is 
to explain how eq. (62) was used in generating Fig. 4. 

If p is the mean occupancy of the server (p = E(x,)) and a is the offered 
load in a GI/M/1-loss system, then it is easy to show that the peakedness 6 
2 satisfies 

I-P 



"(*?)■ 



For GI/M/1-loss system the call congestion is 0(p) and is related to the 
time congestion (p) by a(l — <p(p)) = p. So using z(p) = 1/[1 — 0(p)] — 
a yields the result. 2 ((/>(•) is the L.S. transform of the interarrival time 
distribution.) Therefore, specifying p and z uniquely determines a. 
Hence, with a and z known, we obtain the equivalent random parameters 
and use the three-moment match to obtain the switch parameters (see 
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Ref. 8). Using this procedure, we obtain R(t) vs p parameterized by z. 
Equation (1) then yields Pi,o. 

APPENDIX B 

The Likelihood Statistic Based on the Observable Data 

For ease of derivation, the likelihood statistic derived in Section 4.2 
was based on the raw data x m = (x h - • ;x m ) rather than on the observable 
data [t(m),n{m)\; t{m) and n(m) are defined in Section 2.1. We will now 
study i m (t,n) and verify that the two statistics differ only in their end- 
effect terms. We will also examine the end-effect term based on t(m) and 
n(m) and show how it "tracks" the end-effect term based on x\ and 

We begin by expressing the probability of x m = {x it - • %x m ) in terms 
of t(m), n(m), x\, and x m : 

Lemma 3: If [Xi] is a binary state stationary Markov chain with tran- 
sition probabilities Po.i and Pio and if x m = {x\,- • ;X m ), then 

P(x m ) = P( Xl ) x Q(x lt x m ) x pfJRi - P li0 )"- (£/2) 

XP©l-Po,i) m """ ( ' /2) . (64) 

where 

Q(x h x m ) = Pftr m,/2 Pjr[* 1+ * m)/2 Po,[ xl - Xm)/2 P& + *™> /2+1 (65) 

and 

p if x\ = 1 
1-p if *! = ()' 



P(*i) = 



Proof: This result is obtained from lemma 2 using eqs. (12a) and 
(12b). 

For convenience, we introduce the following notation 

(i) / (-;n,P h0 ) = P$U - Pi,o) n-( ' /2) - Iff is even then of course 

fc(^7l,P li0 ) 

\* I i n \ 

W 



i-n-t/2 



(ii) g (^;m - n,P ,i) = P$(l - Po,i) m - n - 

(Hi) P x<y (t,n) = P(t(m) = t,n(m) = n,x\ = x,x m = y). 
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(iu) S XA .(t,n) = number of binary m-tuples satisfying t(m) = t, n(m) 
= n,x\ = x, and x m = y. 

Thus, we can write 

P(t(m) = t,n(m) = n) = £ £ P x , y (t,n) (66) 

x = 0y=0 

and 

Px^(i»") = f (^ n ' P Lo) # (| m " n ' P °-i) P(*)<S«or('.»)Q(*0'). (67) 

Equations (66) and (67) imply that 
P(t, n(m) = t, n(m) = n) 



= f (!»» P i.o) E (|™ " w,Po.i){ £ LPU)S x , y a,n)Q(x,y)j 
and therefore it is easily seen that: 

, . , . P*(t(m) = t,n(m) = n) 
e m (t,n) = log — — -, 

P(t(m) = t,n(m) = n) 

= \a--an\ + \fi^-b(n-m)\ + E(t,n), 
where 



(68) 



(69) 



£ £ P*(x)S XO p(t f »)Q*(xjf)' 

M/,/n - log J^y^ J. (70) 

£ £ P(x)S Xty (t,n)Q(x,Y) 

x=0y=0 

Comparing eqs. (70) and (22), we see that $ m (t,n) and ${x m ) differ 
only in their respective end-effect terms. The following result is perhaps 
a bit surprising: 

Lemma 4: £(t,n) = S(x m ) if t is odd. 

Proof: todd=>xi^x m ^ e(0,l) = e(l,0) = a + 6/2 (see eq. (22b)). todd 
=> S(io(t,n) = Si t \(t,n) = so E(t,n) for t odd may be written: 

, (1-p)So,,Q*(0,1) + pSi,oQ«(1,0) 
{ ' n) g (l-p)S .,Q(O f l)+pSi.oQ(l,0) " 

Now using (4) E(t,n) for t odd can be manipulated into the following 
form 

E{t,n) = log ! ■ . 

(l-p)S o .iQ(0,l) + pS 1 .oQ(l,0) 
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but Po.i/Pi,o = Po,i/P\.o = p/1 "" P which completes the proof. 
If t is even then x i = x m , 

e(0,0) = b 

and 

e(l,l)»a. 

For even t , eq. (70) can be written as 



E(t,n) = log 



Sq.o j* ■S'1,1 



<£&«.+ £&*» 



Thus, for even £, E(t,n) is a complicated function* of £ and n. Note 
however that E(t,n) ^ and 



It can be shown that 



a if S ofi (t,n) = 
b if S 1 , 1 a,n) = 



S,.,(«,n) n 2 



So,o(*,n) t 

m — n 

2 

if (£,n) is such that 5ii and Sq.o are nonzero. 



APPENDIX C 

The End Effect E\ /( x lm ; x, m + , )} 
The end effect 

E\I(x im ;x im +i)\ = H(xim+i) - H(x im+ i/x im ), 

where 

H(x im+ i/x im ) = (1 - p)H(x im +i/0) + pH(x im +i/l) 

and 

H(x; m 4-i), Hixun+i/O) and H(x fm+ i/l) 

are the binary entropy function 7/(x) evaluated at p, P ()i i and Pi,o, re- 
spectively (fi{x) = —x log x — (1 — x) log (1 — x), ^ x s 1). 

Table I exhibits E\I(xi m ;xi m+i )\ and E\£{x m )\ as a function of the 
trunks occupancy p — for a normal trunk, with a single-hour accumula- 
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Table I — End effect as a function of occupancy 





(t = 


100 


seconds) 




(t = 


200 seconds) 


p 


E\I(Xim+iVCi 


m)\ 


E\t(x m )\ 


B\I( 


fun+l»*i 


«)} 


E\hx m )\ 


0.10 


0.04 




-1.36 




0.01 




-0.22 


0.20 


0.05 




-1.60 




0.01 




-0.20 


0.30 


0.04 




-1.50 




0.01 




-0.15 


0.40 


0.03 




-1.21 




0.01 




-0.09 


0.50 


0.02 




-0.84 




0.00 




-0.04 


0.60 


0.01 




-0.47 




0.00 




-0.01 



tion period (m = 36 and 18 for the 100- and 200-second sampling options, 
respectively). The M/M/1-loss model is used with 1/p = 180 seconds. It 
is clear that the mean end effect is negligible compared to the mean 
statistic update. 



APPENDIX D 

Occupancy Formulae for a Random-Selection Loss System 

Theorem: Consider an N server Markouian loss system with random 
selection of idle servers and 

(i) N - 1 servers with mean service rate p. 
(ii) 1 server with mean service rate rp (r > 0). 
(Hi) Mean arrival rate A. 

Let p r and p* denote the mean occupancy of the servers with mean 
service rates p. and rp, respectively, and let £ denote the blocking (call 
congestion). Also let and 0" denote the carried load per server in an 
N server and N - 1 server Markovian loss system, respectively, (as- 
suming all servers have rate p), given an offered lead a = \/p. That is, 
= a[l - B(N,a)]/Nand 0- = a[l - B(N - l,a)]/(N - \),whereB(- , •) 
is the Erlang blocking formula. Then 

Wp> / 1U - < 71 > 

r - (r - 1)0 

(ii) B(N,a,r) = B{N ' a \, [B(- , ■) is the Erlang (72) 

r - (r - 1)0 



blocking formula]. 

,r - (r - l)[N/(N - l)]<f> 



Proof: 



/ r-(r-l)iV/(iV- 1J0 \ 
<»>♦***( r _ (r _ 1)0 )*■ 

/ r-(r- l)0- \ 

w * = ( r - (r -i)0; - 



(73) 
(74) 
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(i) This part of the theorem follows immediately from eq. (34b) 

1 

Pr = 

rN 

1 + r[l -B(N- l,a)] 

a 

by using the well known Erlang B recursion 6 

= —B(N- l,a) 



l-B(N,a) N 
and the expression for the mean group occupancy 

a[l-B(N,a)] 

N =0 - 

(ii) This part of the theorem follows immediately from eq. (34a) by 
using the above expression for mean-group occupancy. 

(Hi) The lower bound part of eq. (73) follows from eq. (34c) by noting 
that &(N,a,r) s BiN.a)* and using eq. (71) of this theorem. Thus, we 
obtain 

, r</> — + (N - l) Pr * a[l - B(N,a)] = N<J> 
r- (r - 1)0 

which can be arranged to yield the lower bound in eq. (73). The upper 
bound in (73) is an immediate consequence of (74) since 0~ ^ 0. Thus, 
it remains to prove (74). 

(iu) We prove this part as follows: Using eqs. (71) and (72) in the 
conservation eq. (34c), 

rp* r +(N-l) Pr = a(l-B) (75) 

yields 

_ a(r-l)(l-0) + (iV-r)0 
Pr= (JV-l)[r-(r-l)0] * (?6) 

Therefore, eq. (74) holds if and only if 

f t lu-ivj a(r-l)(l-0) + (N-r)0 

[r- (r- 1)0 J X0 = _ i . (77) 

The right-hand side of eq. (77) can be rewritten as: 

rhs = °_ \r(N - I - a) + a - B[N - (a + l)r + a]\ (78) 



t This follows eq. (34a) by noting that a [1 - B(N,a)] £ N. 
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using 

a[l -B(N,a)] 

Similarly, using </> _ = a[l - fi(N - l,a))/(N - 1), the left-hand side of 
eq. (77) can be written as: 

lhs = \r - (r - l)(a/JV - 1)[1 - B(JV - l,a)]\(a/N)[l - B(N,a)]. (79) 

Now using the Erlang B recursion formula in eq. (79) and rearranging 
terms yields eq. (78). Thus, eq. (77) and consequently eq. (74), of the 
theorem is proved. 

APPENDIX E 

The Sign of the Mean Statistic Update 

The mean statistic update of the optimal 5XB group algorithm is given 
by eq. (51a): 

P*(*i) 



E\£(x m )\ = (m - l)|(«Pi i0 - a)p + (b - 0P o ,i)(l - p)\ + log 



P(*l) 



The mean update of the ad hoc and the additional statistics correspond 
to b = $ = and a = a = 0, respectively. The following lemma is needed 
to study the sign of E\£(x m )\: 



Lemma 5: 


://0 


<q 


< p < 1 define 






a 


= i< 


) and a 
-P' 


then 








a 
q<-<p. 



Proof: Consider the log likelihood ratio statistic £(m) for this simple 
hypothesis testing context: the observed process is n(m) = 2 -1] Xj, where 
Xj are i.i.d. Bernoulli random variables with P(xi = 1) given by q and p 
under Hq and Hi, respectively. Therefore, 

2 b(n;m,p) 

t{m) = log — = an — am 

b(n;m,q) 

and hence the mean of £(m) is m(aq —a) under Ho and m (ctp —a) under 
H\. But as we noted in Section 6.2, E\£\ is negative under Ho and positive 
under H\ in a general discrete setting.* Thus, we must have 

aq — a < 

TRUNK-DETECTION ALGORITHMS 975 



and 

ap — a > 0, 

which completes the proof. 

An immediate consequence of this lemma is that 

(i) Pi.o(*,l)< £ 77^<Pi.o(*,r ) 
a(0,roJ 

and 

(u) PoA<l>,ro)<j^<P 0il (<l>,l), 

where a and a are defined by eqs. (41a) and (41b) and b and B are defined 
by eqs. (45a) and (45b). Since r ^ r implies that Pi,o(0,r) ^ Pi,o(0>ro) 
and Po,i(0,/*) ^ Po,i(0.^o). we see that 

sgn(«(0,r o )Pi,o(0,r) - a(0,r o )) = sgn(5(0,r o ) - P o ,i(0,r)B(0,ro)) 

_ [positive if r ^ ro 
1 negative if r = 1. 
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